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I report BESIII preliminary results on: 

1. Measurement of a{e^e~ —)■ DD) at Ecm = 3.773 GeV 

2. Study of the DD production line shape near Ecm = 3.773 GeV 

3. The first observation of singly Gabibbo-suppressed decay, D —)■ un 

4. Measurement of B{Dg —)■ rj'X) and B{Dg —)■ rj'p'^) . 
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1 Hadronic decays of charm mesons 

Studies of hadronic decays of charm mesons play an important role in the under¬ 
standing the weak interactions at the c-sector and provide inputs for the beauty 
physics. Two of samples accumulated by the BESIII detector [1] that are taken at 
Ecm = 3.773 GeV and 4.009 GeV are very useful to study decays of D and mesons. 

The former is the largest annihilation sample in the world to date, 2.92 fb“^ [2], 
that is taken around the nominal mass of '0(3770) resonance which predominantly 
decays into a pair of D mesons. The latter, consisting of 482 pb“^ [3], also produces 
a pair of D^Dg with a sizable production rate {a{e'^e~ —)■ DgDg) ~ 269 pb), 

providing a clean event environment to study decays of 

In this proceeding, I report four preliminary measurements from the BESIII col¬ 
laboration based on the above two e'''e“ annihilation data. The hrst two results are 
studies about D-pair productions at the vicinity of the '0(3770) resonance, a measure¬ 
ment of observed a{e~^e~ —)■ DD) at Ecm = 3.773 GeV and a study of Born-level line 
shape of a{e~^e~ —)■ DD). I then present the first observation of the singly Gabibbo- 
suppressed decays (SGSD), D —)■ un, and end this report with the measurements of 
BiD^ ri'X) and B^D^ ri'p+). 

2 (T(e+e- DD) at Ecm = 3.773 GeV 

Measuring observed a{e'^e~ — DD) allows us to estimate the number of DD pairs 
produced in our sample by using the integrated luminosity of the corresponding 
sample [2]. This can then be used to normalize the measured signal yields to obtain a 
branching fraction. 

As done by the GLEO collaboration [4], we measure the observed cross section 
by a double-tag technique, pioneered by the MARK III Gollaboration [5]. This takes 
advantage of the fact that D-meson production near the 0(3770) resonance is solely 
through DD. 

Reconstructing one D meson in the pair provides a single-tag yield, with a 
hnal state, i. We seek 9 different final states: D^ —(A'“7r+, A'^tt+tt+tt”), 

and D^ —)• (A"“7r’''7r''', A'^tt+tt+tt'’, A'^Tr+vr+Tr", A'’''A'“7r’''). (Un¬ 

less otherwise noted, charge conjugate modes are implied throughout this report.) 
The detail reconstruction criteria can be found in other BESIII publications, such as 
Ref. [6]. 

Ngrp can be written as = Ndp,-B{D —)■ 0-6,, where Njjjj is the number of DD 
produced and e* is the reconstruction efficiency for the decay mode, D ^ i. Similarly, 
one can have ■ B{D —)■ j) ■ ej. When the pair decays explicitly into two 

hnal states, D ^ i and D —)■ j, we have N]jrp = N^p)-B{D —)■ i) -B^D —)■ j) -eij. Here, 
is the double tag yield when we simultaneously reconstruct the two mesons in 
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the final states of i and j. etj is the corresponding reconstrnction efficiency. Solving 
these for Njjjj, one arrives at; 

iV - = 

■ e, • e, ■ 

The observed cross section is readily obtained by dividing Nf)j^ by the total integrated 
Inminosity. 

We obtain Ngrp from distributions of beam-constrained mass, Mbc, defined as 
^Bc = -\Pd?‘- Figure 1 shows fits to Mbc distributions based on singly 

tagged events for the 9 different final states. We use a signal shape predicted by 
Monte Carlo (MC) simulation. Each of these are convoluted with a Gaussian to 
take into account a discrepancy in resolution between data and MC, while using an 
ARGUS background function [7] to represent the background component. 
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Figure 1: Fits to Mbc distributions of singly tagged events based on the entire ■0(3770) 
sample. Red curves represent the overall fitted shapes and blue dashed curves corre¬ 
spond to the fitted ARGUS background functions. 

As for obtaining A’ot) ^ two-dimensional space, Mgc vs M^q. Due to 

the small background of the doubly tagged events, we simply count the yields after 
using the sidebands of Mbc to estimate backgrounds. 

Averaging the resultant observed cross sections over different final states {D —j 
and D ^ j ), we have our preliminary result shown in Table 1. Our cross sections 
are consistent with the ones measured by the CLEO collaboration [4] . We expect our 
final results to be dominated by systematic uncertainties. 
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Experiment a(e~^e —)■ D^D^) (nb) cr(e’''e —)■ D^D ) (nb) 

This work 3.641 ± 0.010 2.844 ± 0.011 

CLEO [4] 3.607 ±0.017 ±0.056 2.882 ± 0.018 ± 0.042 

Table 1: Comparison of the measured cross sections beween the BESIII preliminary 
results and the ones measured by the CLEO collaboration. Only statistical uncer¬ 
tainties are shown in the BESIII results. 


3 Line shape of cr(e+e DD) 

In the previous section, I report our preliminary result of observed cross section, 
(j(e+e“ —)■ DD), at Ecm = 3.773 GeV. It is of great interest to examine this production 
line shape near the nominal mass of '0(3770) resonance. This is done using the BESIII 
scan data which was taken in 2010, along with the main on-resonant 0(3770) sample, 
in a range of 3.642 < Ecm < 3.890 GeV with the total accumulated luminosity of ~ 70 
pb“^. Such a line shape distribution allows one to extract the 0(3770) resonance 
parameters. Table 2 shows some of the recent experimental measurements on the 
nominal mass of 0(3770) resonance. There is a dehnite (and expected) shift in the 
mass when an interference effect is taken into account. 

Experiment Mp( 377 o) (MeV/c^) 

BES (2008) [8] 3772.0 ± 1.9 

Belle (2008) [9] 3776.0 ± 5.0 ± 4.0 

BABAR (2007) [10]t 3778.8 ± 1.9 ± 0.9 

BABAR (2008) [11] 3775.5 ± 2.4 ± 0.5 

KEDR (2012) [12] t 3779.20^:^l[]:f0g| 
t includes interference 

Table 2: Recent experimental measurements on the mass of the 0(3770) resonance. 


To obtain the resonance parameters, we follow the procedure carried out by the 
KEDR collaboration [12] in which we assume that there are two sources that produce 
DD hnal states: one from the decay of 0(3770) and the other from non-0(377O) 
decays. To represent the non-0(377O) decays, we form its amplitude as a linear 
combination of a constant term, which represents the possible contributions from 
higher cc resonant states such as 0(4040), and a Breit-Wigner form, that corresponds 
to the 0(3686) tail above the DD mass threshold [13]. This approach is known as a 
Vector-Dominance Model (VDM), but we also try an exponential form, instead of the 
Breit-Wigner form, to see how much an alternate form affects the resultant 0(3770) 
resonance parameters. 
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The Born-level cross section, abcn-n, and experimentally determined observed cross 
section, aobs, are related as: 

(^obs{W) = J zdd{WVI - x)aborn{WVl -x)Fisr{x, W‘^)dx. 

Here, z^d Si factor for the coulomb interaction for D'^D~, Fjsr{x, W^) is the ISR 
radiator [14], and G{W,W') (a Gaussian) is there to take into account the beam 
spread at the initial Ecm = W. More details can be found in Ref. [12], 

We extract abomiW) based on aobs(}V) with the above relation. aobsiW) is based 
on the singly tagged events by htting to two-dimensional space, AE vs M^c, where 
AE = Ed — Ebeam witli both signal and background shapes are hxed based on MC 
samples. As an example. Fig. 2 shows projections onto the Mbc axes of such two- 
dimensional hts at Ecm ~ 3.7735 GeV (left) and Ecm ~ 3.7984 GeV (right) based on 
the sum of the three decays (see the 3^*^ column of Fig. 1). Notice that the left 
plot of Fig. 2 peaks at nominal mass of D^, while the right plot of Fig. 2 has a 2”'^ 
peak on the higher side. This is due to the larger ISR effect at this particular Ecm, 
which our MG-based signal shape (green) reproduces quite well. 




Figure 2: Projections onto the Mbc axes (in GeV/c^) of the two-dimensional hts {AE 
vs Mbc) at Ecm ~ 3.7735 GeV (left) and Ecm ~ 3.7984 GeV (right) based on the sum 
of the three decay modes. The blue histograms represent the overall hts, while 
dashed red and solid green histograms correspond to the htted background and signal 
shapes, respectively. 

From these hts at each Ecm, we construct the spectrum of the observed cross 
section, aobs- As an example, we show aobs distribution for the case of D^D~ (red 
points) in Fig 3. There, the solid blue curve is the htted shape to aobs, while the 
corresponding a Bom is represented by the dashed brown curve. The dashed orange 
and green curves are the htted resonant and non-resonant components (here, we use 
the VDM to represent the non-resonant component). 

Table 3 shows our preliminary results on the nominal mass, total width, electronic 
partial width of the "0(3770) resonance. The 4*^ column shows x Bdd^ where 

Bdd = H('0(377O) —)■ DD). This is because our ht is only sensitive to the product of 
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Figure 3: Observed cross section, aobs is plotted in the red points based on the D^D~ 
events in which decays into the 6 different hnal states (see the 1®* and the 2”'^ 
columns of Fig. 1). The corresponding aBorn curve is shown in dashed brown. 


Source 

M^(3770) (MeV/c^) 

ry(3770) ^MeV) 

X Boo (eV) 

BESIIIvdm 

3781.5 ±0.3 

25.2 ±0.7 

230 ± 18 

BESIIlExponential 

3783.0 ±0.3 

27.5 ±0.9 

270 ± 24 

KEDR[12] 

PDG[15] 

3779.31};? 

3773. ± 0.3 

oc q+4.4 

^ 0.0 _Q Q 

27.2 ± 1.0 

160tll 4201^2 

[262 ± 18] X Bdd 


(a) Two solutions were obtained from their fit. 


Table 3: BESIII preliminary results based on the two different forms of the non- 
■^(3770) amplitudes, VDM (-^(3686)) and an exponential shape, are shown, along 
with the result from the KEDR collaboration as well as the current PDG value. In 
the column, B^d = -B('0(377O) —)■ DD). 


the two, but not individually. Our preliminary result is consistent with the KEDR 
measurement. In Tab. 3 we also show a result based on the exponential form to 
represent the non-'^(3770) amplitude. As can be seen, this would likely be one of the 
dominant sources of the systematic uncertainty. 

4 D ^ (jJTT 

For Cabibbo-suppressed charm decays, such as the yet to be observed SCSD D —>■ uti, 
measurements are difficult due to low signal statistics and high backgrounds. For 
the case of D —>■ uir, the most recent experimental search was carried out by the 
CLEO collaboration [16]. They set upper limits, B{D^ —)■ a;7r+) < 3.0 x and 
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—)■ cjTT^) < 2.6 X 10“^ at 90% confidence level (C.L.). In the mean time, H. Y. 
Cheng and C. W. Chiang predict the B{D —)■ ujtt) could be at an order of 1 x 10“^ [17]. 

We start with reconstructing one of the DD pairs with the same 9 hnal states (see 
Fig. 1). Then in the other D decay, we look for —>■ where u —>■ 7 r’'' 7 r“ 7 r° 

and 7 r° —)■ 77 . To improve the signal-to-noise ratio, we also select a certain range on 
the helicity-like angle of cu, ^heiicity, which is dehned as an opening angle between the 
direction of the normal to the u —)■ Tr+Tr^vr® plane and the direction of the parent D 
meson in the u rest frame. We require |i7^| = | cos 6'heiicityl > 0.54(0.51) for D+ (D^) 
that are optimized based on a MC study. 

With additional requirements on Mbc and AE to be consistent with a DD pair 
production, we extract our signal yields by htting to the distributions of invariant 
mass of a; —?• 7 r’'' 7 r“ 7 r° as shown in Fig. 4. We use MC-based signal shapes, along with 
polynomials to represent their background shapes. Figure 4 also shows the expected 
peaking backgrounds (represented by hlled histograms) which are estimated by the 
sidebands of Mbc distributions. The extracted signal yields correspond to a statistical 
signihcance of 5.4(j(4.1cr) for D'^{D^) —)■ respectively. 



M3,(GeV/c^) 



Figure 4: Distributions of invariant mass of cu —>■ 7 r’'' 7 r“ 7 r° for D~^ —)■ 7 r+ 7 r“ 7 r° 7 r’'' (left) 
and D^ —)■ 'n'~^7i~7r^7T^. The solid red lines are the overall hts, while the dashed blue 
lines represent the htted polynomials The hlled histograms represent the peaking 
backgrounds, estimated by the sidebands of Mbc distributions. 


We also check to see if the D —)■ urn candidates produce the expected distribution 
of the helicity angle. Figure 5 shows the distributions of \H^\ in which we can see the 
expected = cos^ 6 'heiicity 

In Fig. 4, we can also see peaks that correspond to D ^ rjTi candidates. We extract 
these candidates by htting to the same invariant mass distributions of a; —?• 7 r’'' 7 r“ 7 r° 
with much narrower ht ranges, and without the requirement on the \H^\. Figure 6 
shows such hts from which we also measure B{D ^ rjn). 
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Figure 5: Efficiency-corrected signal yields in the \Hi^\ bins for candidates of —)■ 

a;7r+ (left) and —)■ (right). The black lines are the htted quadratic shapes. 




Figure 6: Fits to distributions of invariant mass of a; —?• for the candidates of 

D~^ —^ T]n~^ (left) and —)• rjTi^ (right). The hlled histograms represent the peaking 

backgrounds which are estimated by the sideband regions of both signal and tag sides 
of Mbc distributions. 


Table 4 shows our preliminary branching fraction measurements. The measured 
B{D —)■ rjTi) are consistent with the known values [15], while B{D —)■ un) are measured 
for the hrst time. 


5 rjX and r}p^ 

The situation of B{Dg —)■ r]'p^) is rather interesting. If we sum the all known exclusive 
rates with p' in Dg decays in the PDG [15], we arrive at (18.6±2.3)%, while B{Dg —)■ 
rj'X) = (11.7 ± 1.7)% [18]. Among the Dg decays that involve rj', the largest single 
exclusive rate is B{Dg —)■ ri'p~^) = (12.5 ±2.2)% [19]. However, a recent measurement 
is about a half of it, B{Dg —>■ t^'tt+tt^) = (5.6 ± 0.5 ± 0.6)% [20] which appears to 
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Decay mode 

This work 

PDG value [15] 

D+ -)■ a;7r+ 
UTT° 

—)■ ri7T~^ 

—)■ 77r° 

(2.74 ±0.58 ±0.17) x lO”"^ 
(1.05 ±0.41 ±0.09) X 10“^ 
(3.13 ±0.22 ±0.19) X 10“^ 
(0.67 ±0.10 ±0.05) X 10-^ 

< 3.4 X 10-^ at 90% G.L. 

< 2.6 X 10-^ at 90% G.L. 
(3.53 ±0.21) X 10-3 
(0.68 X 0.07) X 10-3 


Table 4: Preliminary result on the measured B{D —)■ utt). 


solve the inconsistency mentioned above. B. Bhattacharya and J. L. Rosner come up 
with two predictions, B{D^ T]'p+) = (2.9 ± 0.3)% and (1.89 ± 0.20)% [21], while 
F. S. Yu et al. predict (3.0 ± 0.5)% [22] by factorization methods. 

We can use our sample taken at Ecm = 4.009 GeV to measure these branching 
fractions to confirm the recent measurement. At this energy, the is produced in 
a pair. To measure the inclusive rate, Dg —)■ rj'X, we employ a double-tag technique 
in which we reconstruct its tag side in 9 decay modes shown in Fig. 7. From these 
Mbc distributions, the single-tag yields are readily obtained. 



Mac (Gev/c^) 

Figure 7: Fits to Mbc distributions of the selected 9 different final states of 
decays. The red curves correspond to the total fits, while the blue dashed curves 
represent the fitted background shapes by the ARGUS background functions [7]. 

To obtain the double-tag yields, we reconstruct the 9 hnal states of Dg decays 
and look for the other Dg decays in the hnal states with rj' —)■ 'n'~^7T~r]{^ 77 ) based 
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on the remaining particles. If there is more than one rj' candidate, we choose the one 
that gives the minimum — M^i{PDG)\. We £t to a two-dimensional space, 

vs Mbc, to extract the signal yields, where Mbc is the tag side of the beam- 
constrained mass. Figure 8 shows such hts, projected onto the Mbc axis (left) and 
onto the axis (right). We use MC-based distributions to represent the signal 

shape. As for the background shapes, an ARGUS background function [7] is used on 
the Mbc direction, while the smooth and peaking backgrounds on the axis 

are represented by a polynomial plus double Gaussian shapes. 

From this £t, 68 ± 14 events are observed as signal candidates. This translates 
into B{D^ — rj'X) = ( 8.8 ± 1.8 ± 0.5)% which agrees with the known value [15]. 



Figure 8 : Fit to two-dimensional space, vs Mbc, where Mbc is the tag side 

of the beam-constrained mass. Shown here are the fitted result, projected onto the 
Mbc axis (left) and onto the axis (right). Solid red curves correspond to the 

overall £t, while dashed blue and green curves are fitted peaking backgrounds. 

To measure B{Dg —)■ we simply use the single-tag method by reconstructing 

Dg —)■ where p'^ —?• tt+tt®. We require the reconstructed p' mass to be within 

3cr of the known mass [15], the invariant mass M^+^o be within 0.17 GeV/c^ of the 
known p mass [15], and hnally its XE be consistent with zero. 

The signal yield is extracted by htting to two-dimensional space, Mbc vs cos6t,+ , 
where 6^,+ is the helicity angle of the 7 r+ from the p decay. We expect to see cos^ 6^,+ 
for —)■ p'p'^ ^ while —)■ p'n^'w^ events should be independent of 6^,+ . 

Figure 9 shows projections onto the Mbc axis (left) of such two dimensional £t. 
On the right, a projection onto the cos 6 * 7 r+ axis with an additional requirement of 
(1.960 < Mbc < 1.980) GeV/c^ is shown. Signal shapes are based on MG simulation. 
To represent the background shapes, an ARGUS background function [7] is used on 
the Mbc axis, while a hxed non-Hj background shape is employed on the cos 6 *.n-+ 
axis, estimated from the Mbc sidebands. 
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The fit yields 210 ± 50 and —13 ± 56 events for —)■ r]'p'^ and Dg —)■ 

candidates, respectively. We normalize the rate by Dg —)• mode to obtain 

B{Dg —)■ p'p'^)/B{D~g —)■ K^K~'k^) = 1.04 ± 0.25 ± 0.07. Or with the known 
B^D'g —)■ K^K~'k^) [15], we arrive at B{D^ —)■ p'p'^) = (5.8 ± 1.4 ± 0.4)% which 
conhrms the recent measnrement by the CLEO collaboration [20]. We also set an 
npper limit on the non-resonant decay, B{Dg —)■ p'7i~^7i^) < 5.1% at 90% C.L. 




Fignre 9: Projections onto the Mbc axis (left) and the cos6'.,r+ axis with an additional 
reqnirement of (1.960 < Mbc < 1-980) GeV/c^ (right) of the two-dimensional £t. 


6 Conclusion 

Fonr preliminary resnlts on the hadronic final states in the decays of D and 
mesons based on the two recent BESIII samples are reported. The measnrements 
based on the world’s largest annihilation sample taken at Ecm = 3.773 GeV 
provide statistically snperior resnlts than the previons experimental resnlts, while 
the stndy of decays of Dg based on the sample at Ecm = 4.009 GeV shows the 
very clean event environment at BESIII. It wonid be very exciting to pursne onr Dg 
program as the collaboration plans to take a few fb“^ of e’''e“ annihilation sample at 
Ecm = 4.180 GeV in 2015 — 2016, where the prodnction rate of is mnch higher, 
a(e+e- ^ D*g^D^) ~ 900 pb. 
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